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Abstract 

The  turbulent  portion  of  the  planar  fuel  jet,  for 
isobaric  subsonic  flow,  under  fast  direct  one-step 
irreversible  reaction  between  unpremixed  fuel  and 
oxidant,  is  analyzed.  Mean  spatial  profiles  for 
the  dependent  variables  are  found  numerically  and 
analytically  from  the  governing  nonlinear  partial 
differential  equations,  based  on  an  explicit  eddy 
diffusion,  and  on  a mean  rate  of  reactant  con- 
sumption proportional  to  the  product  of  the  mean 
mass  fraction  of  fuel,  the  mean  mass  fraction  of 
oxidant,  and  the  appropriate  local  characteriza- 
tion of  the  mean  rate  of  strain.  Results  from  the 
model  are  qualitatively,  and,  in  many  cases, 
quantitatively,  compatible  with  experimental  data. 

1.  Introduction 

In  many  aerothermochemical  devices,  the 
rapid  burning  of  initially  unmixed,  but  highly  com- 
bustible, gases  is  preferred  for  control  and 
safety.  The  essential  point  is  that  the  rate  of 
consumption  of  reactants,  and  the  generation  of 
products  and  of  chemical  exothermic ity  are  usual- 
ly controlled  by  unsteady,  inviscid,  inertial, 
large-scale  mixing.  Once  the  large-scale  struc- 
ture is  broken  down  (to  permit  molecularly  con- 
trolled small-scale  mixing  and  chemical  mecha- 
nisms to  proceed),  the  formation  of  product  gas 
occurs  rapidly.  In  fact,  to  keep  such  devices  to 
practical  dimensions,  the  intuitively  faster  mixing 
assoc iated  with  unsteady,  turbulent  diffusion,  as 
opposed  to  laminar  diffusion,  is  preferred.  (Hence, 
if  the  flow  is  not  naturally  unsteady,  transition  is 
often  artificially  induced.)  Thus,  the  combustion 
of  unpremixed  gases  in  turbulent  shear  flows,  or 
(more  succinctly)  turbulent  diffusion  flames,  re- 
mains one  of  the  subjects  in  aerothermochemistry 
and  most  worthy  of  intense  study. 


Of  course,  not  all  turbulent  reacting  flows  fall 
into  the  category  of  turbulent  diffusion  flames,  and, 
in  these  other  cases,  more  general  approaches 
(retaining  the  roles  of  finite-rate  chemical  kinetus 
and  of  molecular-scale  mixing,  as  well  as  that  of 
large-scale  turbulent  diffusion)  are  necessary. 

It  is  intriguing  for  modelers  to  attempt  these  more 
complicated  cases,  but,  from  a practical  point-of- 
view,  turbulent  diffusion  flames  remain  one  of  the 
most  recurrent  phenomena  in  combustion.  If  so 
practically  important  and  relatively  simple  a prob- 
lem cannot  be  feasibly  treated,  then,  the  outlook 
for  the  treatment  of  these  more  difficult  cases 
dims. 

The  present  work  is  also  motivated  by  what  the 
authors  believe  to  be  a realistic  appraisal  of 
(1)  design  needs;  (2)  available  experimental  data 
to  serve  as  boundary /initial  conditions  to  initiate 
a prediction,  and  as  comparisons  to  validate  such 
a prediction;  and  (3)  current  physical  understand- 
ing of  turbulent  transfer.  A relatively  rapid  and 
inexpensive  means  of  explicitly  characterizing  the 
alteration  in  system  performance  as  a function  of 
alteration  in  controllable  parameters  is  often  what 
is  primarily  sought  by  the  designer.  Properties 
of  the  turbulence  are  often  not  of  as  much  practical 
interest  as  behaviors  of  the  dependent  variables, 
themselves.  In  fact,  properties  of  the  turbulence 
are  usually  not  known,  and  must  be  arbitrarily 
adopted  at  initial  and/or  boundary  stations  in  com- 
putations involving  higher-order  closures.  Where 
sensitivity  to  such  assumptions  dissipates,  the 
flow  may  well  be  close  enough  to  local  turbulent 
equilibrium  that  explicit  lower-order  closures  are 
adequate.  Another  problem  with  such  higher- 
order  closures  is  that  the  conceptually  systematic 
developments  become  ultimately  vitiated  in 
practice,  because  of  the  need  to  adopt  relatively 
arbitrary  statements  concerning  the  integral 
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length  scale  of  the  turbulence  and/or  the  functional 
forms  of  the  probability  distribution  functions  for 
the  dependent  variables.  Unless  consistency  in 
the  level  of  approximation  of  a coupled  nonlinear 
system  is  maintained,  professions  of  accuracy 
become  moot.  Thus,  intricate  formulation  neces- 
sarily incurs  additional  expense,  but  does  not 
necessarily  achieve  additional  validity.  In  an  al- 
most inevitably  semi -empirical  subject  as  turbu- 
lent reacting  shear  flows,  the  simplest  formu- 
lation, adequate  in  terms  of  available  data  and 
engineering  design  for  a particular,  yet  commonly 
occurring,  case,  i.  e.  , turbulent  diffusion  flames, 
appears  worth  identifying. 

Of  special  interest  here  is  the  application  of 
modeling  of  unsteady  diffusion  flames  to  diffusion- 
type  continuous  chemical  lasers,  in  which  highly 
dissociated  fluorine  reacts  with  bimolecular 
hydrogen  (or  deuterium)  to  form  an  activated 
product.  For  this  particular  case,  a direct  one- 
step  irreversible  chemical  reaction  is  the  actual 
chemical  mechanism,  not  a simulation  (Kerber, 
Emanuel  and  Whittier*),  Chemical  lasers  operable 
at  higher  cavity  pressure  levels  are  of  interest; 
such  higher  levels  still  leave  three-body  reactions 
negligibly  slow  relative  to  the  two-body  mechanism. 
While  higher  pressure  increases  the  degree  of 
unsteadiness  in  the  flow  field,  the  Reynolds  num- 
bers are  small  enough  that  conventional  concepts 
of  fully -developed  turbulence  do  not  apply  to  the 
pumping.  In  fact,  consideration  of  the  large- 
scale  structure  recently  emphasized  in  transitional 
and  unsteady  shear  flow  is  doubtlessly  relevant 
(Davies  and  Yule^:  Roshko^).  Modifications  of  the 
large-scale  structure  at  higher  Reynolds  numbers 
do  not  appear  pertinent. 


For  such  large-scale  structure,  both  a local 
instantaneous  theory  for  unsteady  diffusion  flames 
(Bush  and  Fendell*;  Carrier,  Fendell  and  Marble^; 

Carrier  and  Fendell®),  and  also  a mean  theory  for 

78 

the  time-averaged  description  (Bush  and  Fendell’  I 
Bush,  Feldman  and  Fendell^)  have  appeared  in  the 
literature.  The  model  is  particularly  suitable  for 
those  cases,  such  as  hydrogen-fluorine  cold-step 
chemistry,  in  which  the  first  Damkohler  number  is 
large,  so  that,  instantaneously,  a thin  flame  arises. 
(Large  first  Damkohler  number  implies  that  the 
reaction  rate  is  large  relative  to  the  flow  rate,  so 
that  the  aerothermochemical  system  is  in,  or  very 
close  to,  chemical  equilibrium.)  In  such  cases,  in 
time-average,  the  mean  rate  of  reactant  con- 
sumption is  taken  to  be  related  to  the  product  of 
the  mean  mass  fraction  of  oxidant,  the  mean  mass 
fraction  of  fuel,  and  the  principal  mean  strain 
rate  (for  the  parabolic  formulation  appropriate  for 
the  thin  shear  layers  of  interest  here).  It  is  argued 
that,  in  general,  the  magnitude  of  the  principal 
mean  strain  rate  indicates  the  level  of  local  chemi- 
cal activity  in  a turbulent  diffusion  flame  in  which 
macroscopic  mixing  is  rate-controlling. 


The  explicit  local  algebraic  (as  opposed  to  field- 
type  differential)  expression  adopted  for  the  mean 


rate  of  chemic  al  reaction  is  compatible  u.  it ti  the 
adoption  of  an  explicit  local  algebraic  (eddy  vis- 
cosity) expression  for  the  mean  rate  of  diffusive 
transport.  It  is  believed  that:  ill  it  is  inconsis- 
tent to  adopt  a second-order  closure  for  the  mean 
kinetics  or  eddy  transport,  and  an  explicit  (or 
first-order)  closure  for  the  other:  it  is  pre- 

mature to  examine  second-order  closure  for 
either  the  mean  kinetics  or  the  mean  diffusive 
transport  until  the  local  explicit  algebraic  re- 
lations, into  which  such  formulations  degenerate 
in  local  turbulent  equilibrium,  are  first  examined 
for  plausibility;  and  (3)  the  additional  omplexity 
of  second-order  closures  for  nondefec t -layer-type 
flow  geometries  (Spalding*®'  **;  Ooms  and  Wicks’"; 
Spalding  *^)  has  not  been  justified  by  proof  of 
sufficient  additional  validity  in  prediction,  and. 
for  defect-layer-type  geometries,  explicit  clo- 
sures are  quite  competitive  with  second-order 
closures  (Reynolds**). 

The  present  model  has  been  studied  in  detail 
for  the  mixing  layer  formed  by  two  parallel 
streams,  a geometry  of  interest  towards  the  end 
of  the  lasing  zone  in  a diffusion-type  continuous 
chemical  laser  cavity  (Bush,  Feldman  and 
Fendell^).  However,  whereas  the  fully -developed 
mixing  layer  is  tractable  for  analysts,  it  is  diffi- 
cult for  experimentalists,  and,  thus,  some  em- 
pirical constants  in  the  theory  remain  unassigned 
for  lack  of  data.  Experimentalists  have  treated 
the  fuel  jet  exhausting  into  an  oxidant-containing 
ambient  (e.  g.  , Hawthorne,  Weddell  and  Hottel 
Kremer**1;  Kent  and  Bilger*';  Bilger  and  Beck*®), 
and,  under  the  assumption  that  values  for  empiri- 
cal factors  are  transferrable  between  the  geo- 
metries, a solution  is  sought  here  for  the  two- 
dimensional  fuel  jet.  The  major  new  feature 
introduced  by  the  fuel  jet,  without  counterpart  in 
the  (downstream)  mixing  layer,  is  the  plane  of 
antisymmetry  for  the  principal  mean  strain  rate 
(henceforth  referred  to  as  a line  of  symmetry, 
since  a transverse  cross-section  only  needs  to  be 
considered).  Because  the  mean  strain  rate  is  odd 
about  the  axis,  clearly,  it  is  the  magnitude  of  the 
strain  rate  that  is  relevant  for  modeling  the  mean 
consumption  rate.  The  question  has  been  raised 
as  to  whether  a model  giving  zero  reaction  rate 
at  the  axis  can  be  used  to  describe  measured 
properties  of  turbulent  jet  jor  wake)  flows  with 
diffusion  flames  (Williams  ),  and  this  point  is 
to  be  resolved  below.  Because  this  work  is  moti- 
vated by  interest  in  pumping  phenomena  in  the 
cavity  of  a diffusion-type  continuous  chemical 
laser,  as  noted  earlier,  this  question  definitely 
warrants  an  answer:  the  upstream  end  of  the 
pumping  zone  has  wake-like  character  that  evolves 
into  mixing-layer-like  character  only  towards  the 
downstream  end  of  the  pumping  zone  in  chemical 
lasers  (Grohs^®).  Capacity  of  the  model  to  treat 
a turbulent  diffusion  flame  for  the  planar-jet  con- 
figuration, with  its  plane  of  antisymmetry  for  the 
principal  mean  strain  rate,  and  to  treat  such  a 
flame  for  the  mixing-layer  configuration,  with 
its  plane  of  symmetry,  suggests  that  the  model 
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can  treat  the  turbulent  diffusion  flame  for  a wake 
evolving  into  a mixing  layer. 


It  is  worthwhile  to  recall  some  well-known 
properties  of  turbulent  jets,  in  order  to  delineate 
the  portion  of  the  flow  field  of  pertinence  here. 
When  fluid  emerges  from  a nozzle  into  the  stagnant 
ambient,  marginally  unstable  mixing  layers  form 
at  the  lips  of  the  nozzle,  and,  eventually,  merge 
on  the  plane  of  symmetry,  where  the  constant- 
velocity,  constant-thermodynamic -property  'po- 
tential core'  terminates  (Fig.  1).  This  distance 
is  on  the  order  of  six  nozzle  widths  downstream 
(Beer  and  Chigier2*;  Davies‘S),  and  methods  for 
calculating  this  distance  have  been  advanced 
(Lessen  and  Paillet^  ).  This  distance  from  the 
exit  is  sometimes  referred  to  as  the  'breakpoint', 
because  the  jet  begins  spreading  at  a greater  rate 
from  this  distance  downstream.  Actually,  even 
after  this  merger  distance,  an  (additional)  dis- 
tance of  several  nozzle  widths  is  required  before 
the  transitional  flow  becomes  fully  turbulent  and, 
in  the  chemically  frozen  case,  a fully-developed 
turbulence  and  self-similar  representations  of 
velocity,  species,  and  temperature  profiles  be- 
come appropriate.  It  is  common,  since,  in  any 
transverse  downstream  plane,  in  the  chemically 
frozen  case,  the  flux  of  momentum,  of  discharged 
fluid,  and  of  discharged  enthalpy  above  ambient 
are  invariant,  to  treat  the  jet  as  if  it  were  emitted 
from  an  ideal  source,  rather  than  from  a nozzle 
of  finite  dimension.  Since  the  interest  here  is  in 
chemically  reacting  flows,  source  ideas  are  modi- 
fied at  convenient.  In  fact,  in  effect,  it  is  taken 
below  that  initial  profiles  for  the  dependent  vari- 
ables are  given  at  the  'breakpoint1.  In  the  ab- 
sence of  data,  plausible  interpolations  between 
nozzle  properties,  which  hold  at  the  axis  at  the 
merger  point,  and  ambient  values,  which  hold  far 
from  the  axis  at  all  downstream  distances,  are  to 
be  adopted. 


There  are  three  observations  about  the  model 
just  presented  that  warrant  discussion.  (1)  It  is 
being  implicitly  assumed  that  the  distance  to  tran- 
sition is  also  approximately  the  distance  to  ini-^ 
tiation  of  burning.  This  model  is  conventional, 
and  jet  diffusion  flames  are  commonly  described 
as  'lifted',  because  they  stand  off  a finite  distance 
from  the  nozzle  exit.  Attempts  to  describe  this 
stand-off  regime  in  detail  inevitably  entail  arbi- 
trary assumptions,  the  validity  of  which  are  diffi- 
cult to  evaluate  (T  amanini^).  (2)  Although  more 
experimental  data  are  available  for  axisymmetric 


Wohl  and  Shipman^*  comment:  "...  there  is 
always  some  distance  ...  above  the  port  for 
which  the  flame  appears  free  of  turbulence. 
Since,  in  that  region,  mixing  probably  occurs  by 
molecular  diffusion,  and  since  molecular  diffu- 
sion is  a slower  mixing  process  than  eddy  diffu- 
sion, it  is  assumed  that  little  combustion  occurs 
in  this  region.  " 


turbulent  jet  diffusion  flames  than  for  planar  turbu- 
lent jet  diffusion  flames,  the  planar  case  is  being 
developed  here  for  two  reasons.  One  is  that  exis- 
ting diffusion-type  continuous  chemical  lasers, 
ultimately  motivating  this  work,  involve  a planar 
geometry.  1 he  other  reason  is  that,  for  simpli- 
city, a self-similar  representation  of  the  flow 
field  is  employed,  in  order  not  to  complicate  the 
energetics,  which  is  of  prime  concern.  Such  a 
self-similar  representation  has  a greater  domain 
of  validity  for  the  planar  jet  than  for  the  axi- 
symmetric jet.  (3)  The  linear  rate  of  planar  jet 
growth  with  downstream  distance,  postulated  by 
self-similar  analysis  for  the  fully -developed  turbu- 
lent flow,  has  been  questioned  recently  as  perhaps 
slightly  in  error  far  downstream  i Kotsovinos1- ' I. 
Even  if  such  a discrepancy  does  exist,  both  the 
magnitude  of  the  effect,  and  the  fait  that  distances 
far  downstream  of  the  flame-tip  length  are  of  less 
interest  make  this  observation  of  little  significance 
for  present  purposes. 


II.  Formulation 


The  direct  one-step  irreversible  bimolecular 
chemical  reaction 

F + O - P,  ( 1 ) 

in  which  fuel  F and  oxidant  O yield  product  P, 
is  considered.  The  flow  geometry  studied  is  that 
of  a low-speed  symmetric  planar  jet  of  one  reac- 
tant (say,  fuel)  exhausting  into  a stagnant  atmo- 
sphere containing  the  other  reactant.  In  the  limit 
of  (1)  13  = P u t /p  , the  (reference)  Reynolds 

number,  gofngF  to  infinity,  and  (2)  M u /(yRT 

r r 
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Gutmark  and  Wygnanski  remark:  I he  nor- 

malized turbulent  intensities  on  the  [c  enter-] 
plane  of  the  jet  attain  their  self-preserving  state 
about  thirty  slot  widths  downstream  f rom  the 
nozzle.  ...  The  approach  to  self-preservation  in 
a two-dimensional  jet,  thus,  occurs  much  earlier 
than  in  an  axisymmetric  jet  ...  or  a two-dimen- 
sional wake.  In  contrast  to  an  axisymmetric  jet, 
there  is  no  indication  here  that  the  transverse 
and  lateral  components  ...  of  the  velocity  fluc  - 
tuations attain  self-preservation  long  after  the 
[streamwise]  component  does.  " 

This  variance  in  the  linear  rate  of  growth  is 
ascribed  to  the  nature  of  the  flow  outside  the  jet 
induced  by  the  jet,  itself. 

*** 

Since  th^exper iments  of  Hawthorne,  Weddell 
and  Hottel,  it  has  been  established  that  the 
ratio  of  flame  length  (by  any  reasonable 
phenomenological  definition)  to  initial  jet 
radius  is  invariant  with  initial  jet  velocity 
in  turbulent  flow,  to  excellent  approximation. 

In  nondimensional  terms,  the  ratio  is  invariant 
with  Reynolds  number,  for  Reynolds  numbers 
(footnote  continued  on  next  page) 
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the  (reference)  Mach  number,  going  to  zero, 
the  (nondimensional)  boundary-layer  approxima- 
tions to  the  governing  conservation  equations 
for  the  time-averaged  description  of  the  turbulent 


'eacting  jet  a re: 

(Pu).  + (Pv), 

0: 

(2) 

X y 

puu,  f Pvu.  = 
x y 

(P<u.  ),  ; 

y y 

(3) 

PuT,  + PvT,  = 

JL  (Pf  T.  ) 

, + Qw: 

(4) 

x y 

ffh  y 

y 

PuY.,  t P v Y . , 

= » (P*  Yj, 

1 

),  + w.. 

ix  i y 

y y t 

with  i = F,  O,  P.  (5) 


Here,  (1)  x = x /■*  and  y = y . / 1 , with  x 2 0, 
k r k r 

0 s y *•  x , are  the  streamwise  and  normal  spatial 
coordinates,  respectively,  with  * = reference 

length  (the  initial  jet  width):  (2)  u = o^/ur  and 
v = v,  / u r are  the  mean  streamwise  and  normal 
velocity  components,  respectively,  with  ur  = 
reference  velocity  (the  initial  jet  velocity); 

(3)  p - Pyl p r and  T = Tk/Tr  are  the  mean 
density  and  (absolute)  temperature  of  the  mixture, 
respectively,  with  PrRTf  = pr  , where  pf  = 
reference  pressure  (the  background  and  jet 
pressure);  (4)  Y.  is  the  mean  (stoichiometrically 
adjusted)  mass  fraction  of  species  1; 

(5)  Q = Qk/c  T is  the  specific  heat  of  com- 
bustion; whift  r w (=  -w^.  r -Wq  = ipl  = 
wk/(Pru  /*  ) is  the  mean  rate  of  reactant  con- 
sumption (and  product  generation);  (6)  < = 
is  the  eddy  viscosity  (or  turbulent  momentum- 
diffusion  coefficient);  while  <7^  and  (7^  are  the 
(constant)  turbulent  Prandtl  number  and  turbulent 
Schmi  numbers  for  species  i,  respectively, 
such  that  f/<7u  and  </(T  are  the  turbulent  heat- 
conduction  and  species-diffusion  coefficients. 
Further,  in  this  formulation,  the  (nondimensional) 
equation  of  state  is  postulated  to  be: 

P = PT  = 1,  (6) 


(footnote  continued  from  previous  page) 

in  excess  of  a couple  thousand.  This  result 
infers  that  molecular  diffusion  is  negligible 
relative  to  inviscid,  inertial,  macroscopic 
processes  for  gaseous  turbulent  jet  diffusion 
flames.  Thus,  molecular  diffusion  is  neglected 
relative  to  eddy  diffusion  in  the  formulation. 

* Here,  the  boundary-layer  approximations  are 
invoked  on  an  intuitive  basis.  That  these  ap- 
proximations are  valid  ones  is  not  at  all  a 
straightforward  matter  to  demonstrate  for 
turbulent  flows.  It  must  be  proved  experi- 
mentally (or  otherwise)  that  any  terms  dropped 
from  the  equations  are,  in  fact,  small  compared 

to  the  terms  retained. 


where  p Pk/Pr  thp  mean  pressure,  and  the 
molecular  weights  of  all  species  present  are  taken 
as  compa  rable. 

In  the  analysis  that  follows,  the  mean  reaction 
rate  is  postulated  to  be: 

w - (P0 1 u,  |)  Y Y (7, 

y F O 

where  d is  the  (nondimensional  I function  that 
characterizes  the  rate  of  chemical  consumption 
relative  to  the  rate  of  species  transport.  Because, 
for  the  case  of  predominant  interest  here,  in 
which  y = 0 is  an  axis  of  (even  or  oddl  symmetry 
for  all  the  dependent  variables,  the  factor  u,  is 
odd  about  y - 0,  the  expression  for  w is  con- 
tinuous at  y = 0 , but  w,  (and  all  higher 
derivatives)  are  discontinuous.  The  source  of 
this  difficulty  is  that,  more  generally,  | u,  | 

is  , where  $ is  the  velocity-gradient  factor 
of  the  conventional  laminar  dissipation  expression 
written  in  terms  of  the  mean  velocity  field  com- 
ponents. Whereas  ; | u,  | is  a satisfactory 
approximation  to  throughout  most  of  the  flow 

field;  near  y = 0 (where  u,  -•  0)  , other  dis- 
carded terms  in  <*9  should  be  restored  in  the 
leading-order  approximation  for  the  local  strain 
rate.  For  simplicity,  a conventional  parabolic 
formulation  is  adopted  here,  and.  near  y 0, 
an  even  function  of  y,  which  recovers  whatever 
degree  of  smoothness  is  believed  important,  may 
be  'patched'  into  the  (elsewhere  uniformly  valid) 
expression  for  | u,  | . The  amount  of  chemical 
activity  at  the  axis  is  believed  to  be  of  such 
magnitude,  relative  to  the  amount  of  chemical 
activity  off  the  axis,  over  most  of  the  axial 
positions  of  interest,  and  the  nature  of  the  experi- 
mental data  for  a reacting  jet  is  such,  that  it  is 
believed  that  adoption  of  (7)  is  currently  accept- 
able on  physical  grounds.  Thus,  (7)  is  modified 
locally  at  y = 0 only  as  required  analytically. 

For  R ® , M_-»  0 , and  pi,  with  the  intro- 
duction of  (1)  the  coordinate  transformation 

y 

(x,  y) -•  (x,  z),  with  z = ^ (y^dy',  <8at 

and  (2)  the  velocity  transformation 

(u,  v)  "•  (U,  V),  with  U = u,  V = u z,  ♦ v z,  , 

x y 

(8b) 

(2)  - (7)  may  be  re-expressed  as: 

U,  f V,  = 0:  U = *,  . V x - *,  ; (9) 

X Z Z X 

♦ 

A preliminary  suggestion  is  that,  in  general,  d 
is  inversely  proportional  to  m *(V  mr.  **Qm0). 
where  m.  and  tA  are  the  molecular  weight  and 
stoichiometric  coefficient,  respectively,  of 
species  i.  Here,  the  case  , P = I,  m^,  m^.s 

m it  being  examined,  so  that  m •»  2m  . 

A X 


4 


UU,  + VU,  - (DU,  ),  = 0; 

X z z z 

UT,  + VT,  - — (DT,  ), 
x z z z 

= Q(B|U.zl)  YfY0; 
UY.,  + VY.,  - JL  (DY.,  ), 

IX  l Z C7  . I Z 

1 


z 


(10) 

(11) 


= -<b|u.zI)YfY0.  with  i = F,  O.  (12) 

In  the  above  equations,  D = t/T^  and  B = |3/T. 
For  this  system  of  equations,  it  is  taken  that 
stagnant  uniform  conditions  characterize  the 
ambient,  specifically: 


From  (11)  and  (12),  it  is  found  that  th.  . onser- 
vation  equations  for  these  Shvab-Zelduvich 
functions  are: 


U4>  , + V<t>  , , I,  0. 

q x q z O q z z 

Directly,  it  is  seen  that  the  following 

holds: 


I 17) 


i nt  eg  ra  I 


fc  = aY  + b and/or 


T + £ Q(Y  + Y ) a(Yr-YJtb,  i 1 8a  I 

F O t*  CJ 

where,  from  (15), 


U -*  U = 0,  T - T = 1, 


yf^yf~=0'  Yo~  Yo»  =con8t- 


(13) 


For  symmetric  initial  conditions,  at  the  center- 
line,  it  is  taken  that 


V - 0,  U,  , T,  , Y_,  , Y , - 0 

z z F z O z 


as  z "*  0 . 


(14) 


(Obviously,  for  the  sake  of  specificity  only,  the 
jet  is  taken  to  contain  fuel,  but  not  oxidant;  and 
the  ambient  is  taken  to  contain  oxidant,  but  not 
fuel.  ) The  initial  conditions  at  the  jet  exit  plane 
are 


U-  U = 1,  T - T = const.  ■*  1. 

Yf~  YFj=  const.,  Y0^Y0.  =0 
for  0 * z < z = £ 

1 j 

as  x — x.  = 0;  (15a) 

U - U_  = 0,  T - T,  = 1. 


yf~  yf-=0-  yo~  yo«  = con,t> 

for  * =*  (y~)<  z<- 

' i 

as  x - x^  = 0.  (1  5b) 

The  specified  functions  of  (15)  are  compatible 
with  (13),  and  with  the  requisite  of  initially 
unmixed  reactants. 


28 

Since,  experimentally  (Batt  ),  0.  , V.  “O,  it 
is  convenient  to  introduce  the  (turbulent)  Shvab- 
Zeldovich  functions  ♦q  , linear  sums  of  the 
dependent  variables  T and  Y^.  The  functions 
adopted  here  are: 


V(Yf-V'Ys 

♦T  * T ♦ * Q(Yf  ♦ Yq>  ■ fc. 


(16) 


IT  1)  *■  k QIY  - Y ) 

_J tl 9Z. 

<Y_.  * Y„  I 

Fj  O” 


1 + 


fV1)YQ~fQYF,V 

(YFj  fY0^> 


i 18b) 


It  is  noted  that  (18)  can  be  rewritten  to  yield  the 
following  expression  for  T as  a fun.  tion  of  Y 

and  Y0: 


(T.-l) 


T = 1 + 


yf  + <yo> 


- Yo» 


(Y 


Fj 


+ Y 


Q Y_.Y 


Fj  O* 


(YFj  + YD< 


0- 


o* 

yf 


Y 


O 


Physically, 


yfs 


Fj 
and  Y, 


O 
s Y 


O’ 


(191 


Thus, 


Fj 

if  the  (so-called)  adiabatic  flame  temperature 
(defined  to  be  the  temperature  achieved  where 

both  Y_  and  Y_  are  zero)  is  denoted  by  T ,, 
OF  af 

then  T is  the  maximum  temperature  achievable 
in  the  ffow  system.  Since,  experimentally 
(Hawthorne,  Weddell  and  Hottel  ; Kent  and 
Bilger^),  it  is  known  that  the  turbulent  diffusion 


flame  has  finite  mean  structure,  i.  e. 


O 


and 


Y_  never  vanish  simultaneously,  it  follows  that 

r 


1 * T<  T , = b = 1 

af 


(T  j - 1 ) Y 0lB  ♦ QY  f,  Yq, 


<YFj 


+ Y _ ) 
O® 


(20) 


It  is  noted  that,  even  with  the  obtaining  of  the 
above  integral  relating  the  Shvab-Zeldovich 
functions,  the  dynamics  ((9)  and  (10))  and  ener- 
getics ((11)  and  (12))  remain  coupled.  In  general, 
D and  B,  the  turbulent  diffusion  and  reaction 


If  an  effective  flame  position  for  the  turbulent 
flow  is  defined  by  the  locus  Yq  - Y 5 Xf. . then, 
the  flame  temperature  T^  is  given  by 

T,  = T - QX  : (T  - T ) = QX  > 0 , 
f af  f af  f f 

where  T^.  and  X^  are,  in  general,  not 
constants. 
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coefficients,  respectively,  are  functions  of  the 
temperature  (and/or  mass  fractions  of  the  reac- 
tants) and  the  properties  of  the  turbulence  (in 
second-order  closure).  As  remarked  in  the 
Introduction,  consistent  with  (1)  the  preliminary 
character  of  this  investigation,  (2)  the  nature  of 
the  formulated  model  of  the  mean  rate  of  reactant 
consumption,  and  (3)  the  still  imperfect  state  of 
second-order  closure  for  turbulent  jet  flows,  in 
what  follows,  it  is  taken  that  D and  B are,  at 
most,  specified  spatial  functions. 

111.  Approximation  for  the  Mean  Velocity  Field 


the  velocity  components  and  the  vortu  ity  m the 
domain  0 4 tt<  *,  1 5 ? < x)  can  be  written  as 

U 'i,  = f'(»?), 

Z l 

V -*,  - u • T7K  ' (17  > -i  K.TJi  • : (24al 

x y i 

u,  - 'l',  - (u . / z . ) 5 v if) i 

Z Z Z II 

= y(u./x.)4"  L fV/I  . (24b, 

1 1 

Thus,  the  boundary  conditions  of  (2lbi  can  be 
written  as 


With  D a specified  spatial  function,  the  be-  F,(<”)  = 0,  F(0),  F(0)  = 0:  (25al 

havior  of  the  mean  velocity  field  for  the  jet  geo- 
metry is  governed  by  while,  by  definition. 


u,  + 

X 

< 

N 

II 

o 

U = *,  , V^-*,  , 

Z X 

F'(0)  = 1, 

(25b) 

u u, 

X 

+ vu,  - 

z 

(DU,  ),  =0 

z Z 

In  turn,  the  conservation  condition. 

122),  for 

(0  * z < ®,  0 * x < *); 

(2  la) 

momentum  takes 

the  form 

U - 0 

as  z - ® 

, V,  U,  "•  0 as  z - 0; 
z 

(21b) 

" ■ V (f ) 

J 

U - 1 

for  0 s 

z<zj  ’ 

= 2(„2z.,  f’ 
1 1 

f l2 

|f'(7?)1-  dt? 

U - 0 

for  z.  * 
J 

*6r)<z<“ 

j 

as  x ■*  0. 

(21c) 

■Hi'S 

r i2 

|f'(T))}  dTJ  . 

(26) 

Integration  of  the  momentum  equation  of  (21a) 
with  respect  to  z,  subject  to  the  boundary  con- 
ditions of  (21b)  and  (21c),  yields  the  momentum 
flux  K,  an  integral  flow  quantity,  defined  by 


With  thje  introduction  of  the  'effective1  jet 
exit  plane  x;  (>  0)  , the 'normal  length' function 

z . (x)  = z.(x/x.)  , with  z.  = const.,  and  the 'center- 
1 I l •• 

line  velocity1  function  U + (x)  = u^x/x^  with 

Uj  = const.  , and  with  the  introduction  of  (1)  the 

(so-called)  similarity  coordinates 


i 


z 

**<*> 


■W(i)  ■»(!)• 


(23a) 


where  T = (x^/Zj)  = (constant)  spreading  parameter 
(to  be  determined  empirically),  and  (2)  the  stream 
function 


In  what  follows,  an  explicit  eddy  viscosity 
formulation  for  D is  adopted,  namely: 

D(x,  z)  = D(0  = *■ 

= * (u.z.  )£^  = — (u.x.  )£  . (27) 

11  y ii 

where  K = (constant)  eddy-viscosity  parameter 
(to  be  determined  empirically).  This  form  is  the 
one  most  frequently  used  in  the  analysis  of  two- 
dimensional  jets,  because,  in  general,  it  has 
yielded  reasonable  results  (Schlichting2^).  For 
K and  y related  by  * = l/(4y),  the  boundary- 
value  problem  of  (21a)  - (21b)  reduces  to  the 
self-similar  form 


F (TJ)  +2 


F(T))F*(T))+  (f'<»))} 


F'(-)  = 0,  F(0),  F*(0)  = 0, 


0 (0ST)<»); 
F'(0)  = 1. 


(28a) 

(28b) 


*(x.  z)  * *<C.n>  = U„(C)  zji)  F(r?) 

= (u.Zj)  {*F(r» 

= y (Vil?*F(tl)'  ,23b) 

*The  plane  x = x (and/or  £ = l)  is  also  referred 
to  as  the  'brealipoint'2^  and  as  the  'flame  lift-off 
distance',  2* 


Successive  integrations  of  (28a)  yield 


F'(t))  4 2F(rj)  F'(Tj),  = const.  = 

0; 

(29a) 

f i2 

f t2 

F * (TJ)  + |F(t))r  = const,  = F (0)  = 

(f(»)|  = i. 

(29b) 

Further,  (29b),  itself,  can  be  integrated  to  give 
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F(r))  = tanh  t): 


w ith 


F*  t»?)  = 1 - tanh1"  d, 

F'(T))  ' - 2 tanh  »7  ( 1 - tanh^T))  . (30) 

Based  on  the  solution  of  (30),  it  is  determined 
that 


(r) 


for  z 0 as  x * 0, 


8 'I 


u “ 4(x  / «pi-v3;sv— / •' 


r 1 \ 4 2 4 2 * 

K = 2z.  = (—  ) = — (u.  z.)  =—  (u.  x.)  . (31) 

l \T . / 3 ii  ly  ii 

J 

Thus,  the  solutions  for  'i  and  U,  V are 

'ii  - ^ (u.x.)  4 tanh  t) 
y i l 

/ 3 N 't 

«=  ^ ' j 4 tanh  T)  for  u^  « 1;  (32a) 

j 

U = u.4  ^(1  - tanh2  t)) 

l 

* 4 ^ ( 1 - tanh2  f))  for  u^  «»  1 . 

V = ~-  u.  4 -I  T|(  1 — tanh2  t?)  - £ tanh  T)1 

> i I •) 

**  ^ 4 ^ jf)(  1 - tanh2)?)  . £ tanh  rj  J- 

for  u.  * 1 . (32b) 


The  above  solutions  to  (21a)  - (?'b)  are  not  the 
most  general  ones;  they  are,  in  fact,  the  self- 
similar solutions  for  a jet  issuing  from  a slit  at 
4 = (x/xj)  = 0,  and,  thus,  do  not  satisfy  the 
starting  conditions  of  (21c).  In  light  of  the 
improvement  to  be  expected  from  alternative, 
less  tractable  procedures,  it  is  taken  that 

U « C'^(l  - tanhZt)) 


& (1-unhZ  {(V)(;r)  z}j- 


(33a) 


* .4 

If  it  is  taken  that  U(4.0)  = UJ$)  = u.4  8 for 

1 < 4 < * and  that  U(4, 0)  = 
it  is  consistent  to  take 


for  <r<  4 < 1. 


0)  = 1 for  (T< 

ul“(‘Vl/ur) 


1. 


For  u.  = 1,  from  (31),  it  follows  that  z.  = 
(x./»  = (3/4  Tj)  and/or 
= [(3/4)/((Tk)./Tr)]  y.  Thus,  the  (non- 
dimensional)  location  of  the  'effective'  jet 
exit  (or  'breakpoint')  plane,  x^ , is  linearly 
proportional  to  the  jet  spreading  parameter, 

y . 


for  z > 0 as  x * 0,  (3  3b) 

encompasses  the  initial  conditions  of  interest  for 
present  purposes  ( i . e.  , for  x x . and/or  4 ^ If. 
Thus,  (32)  furnishes  the  solution  for  the  jet  flow 
field,  with  Van  (empirically  I assignable  para- 
meter. 


IV.  Rate  of  lleactant  Consumption 

With  <T  “(T,  the  behavio^r  of  the  (modified) 
Shvab-Zeldovich  function  V (Y  + Y(J  i 

= * Yq>  j for  the  let  geometry  is 

governed  by 


U Y, 

fVY,  - 

- (DY.  ).  0 

X 

z 

<7  z z 

(0  4 z 4 ».  0 4 x ' »)  ; 

134a) 

Y - 0 

as  z "*  * 

, Y,  * 0 as  z - 0; 
z 

( 34b  1 

Y - (Y 

Fj  + Y0»* 

for  0 4 z 4 z - * l— . 

h (~~j  4 z < *>  as  x - 0. 

j 

Y - 0 

for  z.  = 
J 

(34c) 

The  species  flux  integral  associated  with  this 
function  is 


N=2^UYdz  = 2zj(YFj4Y0J 

• (7-)°W-  <35> 

j 

In  light  of  the  approximations  already  adopted  for 
the  velocity  field  (cf.  (32)),  it  is  taken  that  Y is 
of  the  form 

Y(x,  z)  = Y (|,  *?)  = Y.(C)  G(tj)  --C.C^G(T))  . (36) 

v l 

With  introduction  of  this  form  for  Y and  the 
above-mentioned  forms  for  ♦ and  its  derivatives, 
the  boundary -value  problem  of  (34a) -(34b)  reduces 
to 


G'(TJ)  + 20 

("f(t))  <u'(V)  + F#(*l)  G(r?)  ] o 

(0  4 »?  < ®); 

(37a) 

G<»)  = 0, 

G'(0)  = 0; 

(37b) 

G(0)  = 1. 

(37c) 
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Integration  of  (37a)  yields 
g'(i?)  + 2aF(TJ)  G(»?)  = const.  - 0:  (38a) 

n 

G(rj)  = exp  \ -2 a \ F(rj  ) dt?  } 

0 

2 |a 

- l-{F(rj)\  (I  - tanh3rj)Cr.  (38b) 

Based  on  this  solution  for  G,  it  follows  that  the 
species  flux  integral  N is  given  by 


N = 2z . ( Y + Y 

) Fj  O' 


»)  = 0r),YFj+Yo»' 

j 


- 2(u.C.z.)  \ F'(TJ)  G(T))  dV 
1 11  0 
2 00 

= -(u.C.x.)  \ F (11)  C(r))dV 

■ yrKrKf  5 <»> 


Ihis  form  is  not  compatible  with  the  conditions  at 
£ - (x/x.)  0.  As  with  the  velocity  field,  it  is  a 

self-similar  form  compatible  with  a passive 
scalar  issuing  from  the  jet  slit  at  £ 'x/x.l  0, 
Again,  in  view  of  the  improvement  to  ie  expet  ted 
from  alternative,  less  tractable  procedures,  it 
is  taken  that 


V - <YFj+Y0»> 


x c5  / 2l  'Hl  i x •'  / 

-J  ( I -tanh  * 


encompasses  the  initial  conditions  of  interest  for 
present  purposes  (i.e.,  for  x*=  x and  £ ll. 

Consider  now  the  behavior  of  Y . With 
CT.  “ a,  and  with  B a specified  constant,  since 
Y^  = - (Y  - Yp-).  the  boundary-value  problem 

for  Y (cf.  (12)  - (15)  ) can  be  written  as 

F 


ci  = ui(7r)(YFj  + Yo= 


UY  , tVY  , --  (DY  , ). 

F x F z CT  F z z 


~ (irJ  <yfj  * Yo»’  for  “i-1,  (40a) 

where 

00 

rt  2 1 

1 = 1(0)  = \ (1  - tanh  T))  d*) 

J0 

r 2 a 

= \ (1  -IP  ) d<0.  (40b) 

J0 

Values  of  1(a)  are:  1(0)  = 1.  l(fc ) = ir/4, 

1(1)  = 2/3,  1(2)  = 8/15,  ...  , The  range  of  a of 
interest  is  (approximately)  £ *a  * 6/5  (cf. 
Reynolds^).  It  is  noted  that,  for  £ s a * 1 , 
(8/3w ) s (2/31)  s 1 . In  what  follows,  it  is 
taken  that  (2/31)  <»  1 . 

Thus,  based  upon  the  preceding  analysis,  the 
solution  for  Y is  of  the  form 

Y =Ui  in)  <YFjhYoJr*(!-tanhV 

*•  <YFj  + Y Qoo ) 4 * ^ (1  -tanh2*))'7, 

i.e.,  C.  - (Y_.  + Y ) , 
i Fj  0* 


for  «,.  («■)"  ‘* 


If  an  effective  flame  for  the  turbulent _flow  is 
defined  by  the  locus  Yq®Yf  , then,  Y = Yqs 
and 


.4  2 a O® 

$f  (1  - tanh  0f)  « (y  +Y  ) 

Fj  O" 

(footnote  continued  in  next  column) 


-Blu.z1  Yf{yo„-(Y-Yf») 


(0  s z < ®,  0 

S x < »); 

(43a) 

Y "•  0 as  z — 

®,  Y , - 0 

as  z — 0: 

(43bl 

F 

F z 

Y_  - Y_.  for 

0 s z < z . * 

F Fj 

J 1 

J./’ 

Y_  — 0 for  z.-i(— -)<z<®  as  x-*0.(43ct 
F j \Tj  / 

For  Yf  of  the  form 

Yf(x,  z)  = X($,Tj),  (44) 


(footnote  continued  from  previous  column) 
gives  the  shape  of  the  flame.  Thus,  in  the 
present  model,  this  shape  is  explicitly  a function 
°f  Y , Y , and  a only,  and  is  not  a function 

of  T.  or  T^f®  1).  (It  should  be  recalled  that 
all  species  have  been  taken  to  be  of  comparable 
molecular  we ight.  ) Furthermore,  the  length 
of  the  flame  (with  for  *s 

given  by 

1°..  /--EL*  Yo->)2 

f : Y _ J * 

O°o 

It  is  seen  that  the  (mean)  position  of  the  effective 
flame,  as  defined  here,  is  independent  of  the 
model  adopted  for  the  mean  reaction  rate  w , 

The  thickness1  of  this  effective  flame  does 
depend  on  the  model  adopted  for  w . 
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from  (43),  (44)  and  (32),  (41),  the  bounda ry -value 
problem  for  X(£,  t?)  in  the  domain  (0  s f?  < ®, 

1 S £ < can  be  written  as 

X,  m + 2a  <p  x,  n - 4a  (1  -o2)  4 X, 


= (8CT>B)  -<0(1  - IP2)}- 


-Fc 


Yo--'Ci?'*(1 


w th  IP  = (P(T?)  = tanhH:  (45a) 


X-0  as  T)  - ®,  X,  - 0 as  T)  - 0;  (45b) 

X - X^  = fnc(tj)  as  £ — 0.  (45c) 

To  complete  the  formulation  of  this  boundary-value 
problem  for  X(£,  T|)  , the  profile  at  the  'effective' 
jet  exit  plane  (£  = 1),  i.  e.  , X^  = X^V),  must  be 
specified.  Here,  this  initial  profile  is  taken  to  be 

X.  = X.(TJ)  - A.  (1  - 02)a=  A.  (1  -tanh2?))*7,  (46) 

tit  l 


with  A.  = const,  (to  be  specified).  For  this 

initial  profile  for  X,  it  follows  that  the  initial 

profiles  for  Y and  Y^,  i.e.,  Y^A*!)  and 

Y_.(TJ)  , are 

Oi 


Y - X.  - A.  ( 1 - <02)CT; 
Fi  l i 


(47a) 


Y = Y 
Oi  O’ 


(Y.  - Y ) 
i F l 


= Y - (Y.  - X.)  = Y - (C . - A.)  (1  - ip2)*7. 
O®  i i O®  i i 

(47b) 


On  the  centerline  at  the  'effective'  jet  exit, 

Y _ . ( 0 ) = Y°.  and  Y_.(0)  = Y°  . , with 
Fi  Fi  Oi  Oi 

Y°  . = A.,  Y°  = Y_  - (C.  - A.) . (48) 

Fi  i Oi  O®  t i 


o o 

For  Y„«Y„  and  Y„.«*»0,  consistent  w ith  no 
Fi.  Fi  Oi  24 

combustion  upstream  of  the  'breakpoint1,  it 

follows  that  A.  ••  Y and  that  C."t  (Y  t-Y-^) 
(cf.  (41)).  Hence,  (4^5c ) is  taken  \o  be  ^ 

2 <7 

X - X.  «»  Y_.  (1  -ip  ) . 
i Fj 

with  <P  =IP(T))  = tanh  T),  as  £ -*  1 . (45c) 

Here,  it  should  be  recalled  that,  in  the 
development  of  the  present  model  boundary-value 
problem,  (45),  there  are  two  empirical  factors: 

"Y,  the  growth  rate:  and  B,  the  effective 
Damkohler  number.  Thus,  the  model  retains 
the  one  empirical  constant  usually  present  in 
chemically  frozen  flow,  y,  and  adds  a second  one 
for  the  chemistry,  B.  More  specifically,  in  this 
model,  it  is  the  product  XB  that  is  required  to 
characterize  the  reaction  contribution.  (With 


respect  to  the  fluid  mechanics,  it  is  noted  thal 
the  use  of,  at  least,  one  empirical  constant  to 
assign  the  growth  rate  is  implicitly  present  in 
second-order  closures,  as  well,  i 


With  the  determination  of  Y ^ X from  i45|, 
the  solutions  for  Y and  T are  determined 
f rom 


O 


f -*  2 , 

Y_  M -£  < 1-0  » 1 

O® 

f . -4  2 a 

- 1 Y £ s (1  -ip  I - X 
F.t 


[Y0®-(YFjfY0®'r'(1-<fl‘,'l'T 


X : i 49 1 


T *>  r 1 + (T  - 1)£  ( 1 - C2)a  ] 


+ Cl[  Yj_  £*  ( 1 -ip 


2,a-x 


l + f(T  - 1 ) + QY__ . k£  "'  ( 1 - o'1' ) J 

l J Fj  • 


- ax.  i so  * 


Here,  for  completeness,  it  is  noted  that,  as 

t-  i , 


*=  Y (1  - ip2'0  ; 
Fi  Fj 

(51a) 

Oi-  Yo®{1-(1'^'a  = 

(51b) 

« 1 +(T.  - 1)  (1  - ip2|a  . 
J 

(51c) 

V.  Approximate  Solutions  - I 

• the  transformation 

(0,0),  with  0 4 O tanh  *7 , 

(52) 

the  boundary-value  problem  for  X Xl9,(P),  in 
the  finite  domain  (0  s 9 s 1,  0 4 1)  can  be 

rewritten  as: 


(i  -tp  >x.  .-2(i-o)  ipx,  = - 2(7 n : 

(pip  (P 


n»8x, 


e 

(4yB)(p 


'o®('  '{e<1  ■<°2,<Ti)  j 

|yFj{0(1  -<p2)tT1--x 


X-  0 as  IP-  1.  X,0-O  as  iP-0; 

2 O q , 

X - X.  m Y (1  -IP  ) as  9 - 1 . 
t Fj 


:(53a) 

(53b) 

(53c) 
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The  first  and  second  integrals  of  (53a),  taking 
into  account  the  boundary  conditions  of  (53b),  can 
be  w ritten  as 


X, 


in 

l 1 - CT 
(1-0  ) 


o n ( 
■'o  (i-o,)a 


(54a) 


°i  n do, 

S — To 

Jo  ([-op 


d<P2 

0-o2)Ua 


(54b) 


The  integral  1(a)  has  been  introduced  previously 
the  integral  T (GT ) , with  J(I)  7/15,  l|j  | 

= (3w-4)/9 is  introduced  here  for  the 

fi  rst  time. 

The  solution  of  ( 5b)  is  determined  (o  be 


7.(9)  - 


9 exp  ' \(  1 - 9 ) | 


c 

1 + X \ 

9 exp  \(  1 - 9. 1 \ d9 

1 "9  ‘ 1 

|59l 


The  value  of  X at  the  centerline,  X(6,  0)  = X (9),  In  turn,  with  the  evaluation  of  7(9 1 , the  first- 
is,  thus,  approximation  for  X(9,0)  is  found  to  be 


X°  in  \ 


<0 

r2 

n 

dO  j 

do2 

( 54c ) 

^0 

(i- 

2 a 

, 2 1 -a  ■ 

-o,) 

< 1 -<£>  ) 

0 


To  gain  approximate  solutions,  an  iteration 
scheme  for  the  determination  of  the  solution(s) 
for  X(9,0)  from  the  integral  form(s)  of  the 
boundary-value  problem,  (54)  , is  proposed.  To 
initiate  th:  sc  =me,  it  is  taken  that  the  zeroth- 

approxima  ion  r X(0,O)  is  of  the  form 

X(Q)(9fO)  =X°(9)(1-«52)CT 


YFj7(9)  {9(l-«)2)a}  , 


(55) 


with  7(9)  a function  to  be  determined  (subject 
to  the  condition  that  7(1)  = 1)  (cf.  Bush, 
Feldman  and  Fendell^). 


With  the  introduction  of  this  approximation, 
evaluation  of  the  integrals  of  (54c)  yields  the 
following  (ordina  ry)  differential  equation  for 

7(9): 


9 z'(0)  - [(O  - X0)  + X0  Z(0>]  Z<9)  =0; 

Z(l)  = 1 , 

Here, 

a «*  (2?B>  I Y, 

X 

X «.  (2yB)  J Y, 


0‘ 


(2XB) 


(}) 


O®’ 


(2yB)  J (YFj  + \QJ. 


Fj  ’ 


where 


1 


2 a 


I = I((7)  = V (1-oT  do, 

' 0 


J = J(a) 


2a  r 

1 +a  .) 


, [l-(.-«P2»l+a] 

0 ,i-oV*° 


dO. 


(56a) 

(56b) 

(57a) 

(57b) 

(57c) 

(58a) 

(58b) 


X(l)<9,0>  ^ X(0|(9.O)  ri  -B(1)(0,Ot  , (601 

where  Xjqj(9,0)  is  given  in  (55),  and  where 

R(1)(0,O)  =aAj(o)  - '"x  - X 7(9)  ^ 9 A ((o) . (61) 

For  a = 1 , 


A.(O)  = 


Iv'  " 1 +0  ’ 

2 

O 


AJW)  = 1+0-7 


3 3 

- T <P  • 


(62a  i 


(62b) 


while,  for  a = £ , 


A,(o>=- 


O + 


(w/2)  1- ( 1 - O2)^  v - sin  'o 

(1  - oV 


(6  3a  I 


AJ(<0'  3tt-4 


o(i-o2)* 


(tr/2)  { 1 - ( 1-0  - s in  ^ o 


+ 3 


(1  -o2)S 


+ 2 


(d-o2)^. 


(1-0) V 


2 t 

(1  -or 


. (63b) 


It  is  noted  that  A (0),  A (0)  0,  such  that 

R |fl,0)  = R°  (9)  = Further,  at  the  'break- 
point1 plane,  9=1, 

R(l)(l,0)  = R(1)i(0)  = aAj(o)  - (X  - X)A  j(O) 

r 1 

<w(2yB)Y„  1 1 A.(O)  - .'A  (O) 

O®  I J -j 


i af  A jlot 


(64) 


10 


exp  X 


Thus,  R (1,  0)  = R (0)  = 0,  while,  R (1,  1)  = 
wiln  f,  a (weak)  fund  ion  of  a, 
whose  numerical  value  is  less  than  unity.  In 
general,  0 * R ^ . «D)  s f<*  , for  0 4 <0  5 1 . 

For  Y = Y,J0,«?>  given  by 
r r 

Y «=  X 
F (1) 


then,  Yq  = Yq{8,<9)  is  given  by 

Yo“  Yo.-(Y-X(i)) 

- Yo,|  1 - {e'i  -<o2)a}l 


-YFj  t1-Z,1-R,l,) 


J 

} {eu  Vf}  . 


Z(0)~  [l  + O(0)J  , 

where 


u 1 v 

1 + X B ^ exp  ; \i  1 - v dB  j 

Thnt  f the  downstream  behaviors  for 
7 , along  the  centerline,  are 

o n 1 IO 

yf~  lo9  +---  • 


'o  ~ ' O®  " n Fj  + ' O*  * 9 + 


t - i t(c.-i)  (QY,.  1 e 

*•  J Fj 


To  the  same  order  of  approximation,  T = T(0,<0)  is 
given  by 

(T.-l)tQY 

T <*»  1 + — ^ 1 Y - QX 

,YFjfYoJ  <»> 

•»  1 + (T.  - 1)  (9(1  -to2)0}  J 
+ qyf.  {l  -Z(l  -R(1))}{0(1  -<p2)<7} 

“ I+{(T.  - 1)  +QYFj}  {0(1  -«52)a} 

- QYf.  {z(l -R(1))}{0(1  -t©2)*}  , (66) 

For  0 = 1 (i.e.,  £ = 1),  from  (65)  and  (66),  it  is 
determined  that 

YFi..YFj(l-02)a(l-R(,)i);  (67a) 

YorYo.{'-<»-V} 

- Yf.(1-«,2,'7  R(1)l;  ,67b) 


1-  1 f (T.  - 1)  (I  -<p  ) 
‘ J 


♦ QYFj(!-0)  R(1)l.  (67c) 

A comparison  of  (67)  and  (51)  indicates  that  the 
original  initial  conditions  are  not  satisfied  exactly 
by  the  first  iteration  solutions. 

F rom  (59),  it  is  seen  that,  as  0-*O  ((-*), 

Z,0|~  L.00’  f J ♦ 0(0)  ] , (68a) 


where,  for  explicitness,  it  is  taken  that  n - )|  » 

I j 

QYF  j-  is  a positive  quantity.  Indeed,  based  upon 
(68),  it  is  found  that  the  downstream  behaviors  for 
Yq  and  T , in  general,  are 


'r,  ii,-..  ti  p i 

O O00  r j Ox 

T ~ 1 + {(T.  - 1)  +QY„.  d 0(1  -<02)CT  + , ..  . (721 

With  respect  to  the  temperature,  it  is  seen  that, 
downstream,  for  a fixed  value  of  0 . 1 is  a 

maximum  at  (0  = 0 and  decreases  monotom.  ally 
as  <0  increases.  Since  reaction  consumes  Y , 
it  follows  that  Y 4 Yj_.  0(1  -C2!**-  , the  'frozen 
solution(for  w =0),  and1,  in  turn, that  T is  aug- 
mented (or  left  the  same)  when  reaction  occurs, 
relative  to  the  nonreacting  case.  (1  his  result 
might  have  been  anticipated  from  (501;  here,  it 
follows  directly,  ) 

VI,  Approximate  Solutions  - !1 

A more  explicit,  but,  probably,  less  accurate, 
approximate  solution  of  the  nonlinear  boundary- 
value  of  (45)  is  obtained  now  through  the  use  of 

Oseen  linea  r ization  ( Lew  is  and  Carrier^; 

32  9 

Carrier  ; Bush,  Feldman  and  Fendell  ). 

Here,  (45a)  is  approximated  by 

%r,+  2aT'5%-4a«x-c 

*=  (8cryB)(aY0s[>T)2)  X,  <73a> 

with  the  constant  a considered  to  be  an  assignable 
parameter  (that  permits  the  solution  of  this 
equation  to  more  adequately  approximate  that  of 
(45a)).  The  boundary  conditions  are  those  of  (45bl, 
i.  e. , 


In  Oseen  linearization,  terms  are  approximated 
in  a manner  to  retain  their  essential  roles,  but, 
also,  in  a manner  to  generate  analytic  tract- 
ability.  Justification  lies  entirely  in  the 
a posteriori  demonstration  that,  to  an  adequate 
degree  of  accuracy  for  the  purposes  in  mind, 
results  satisfy  the  original  boundary -value 
problem. 


X -•  0 as  q - “,  X,  -0  as  rj  - 0 . (73b) 

V 

In  lieu  of  (45c)  , it  is  required  that 


0 

where 


‘ Yn 

as  £ - 

1 , q - 0; 

X riq 

- Y H 

as  £-  1, 

(73c) 

1 


i 2 a d«5 

H - H(cr)  = \ (1-tanh  q)  d q , 7T~7, 

Jo  ’od-to  ) 


Clearly,  H((T)  is  defined  for  (T  > 0 only. 
W2,  H(  1 ) = 1,  H( 3 / 2 ) = ir/4.  H( 2)  = 2/3  , 


(74) 
H<$)  = 


Under  the  transformations 

* » ^ 

£ * log  it  V - (2<T ) *) . 

x<c . »?>  * x*(C *.»>*). 


(73a)  takes  the  form 

• * * * * **  * 
%*»,*  fn  x*n* ‘2x't*'Q  n x’ 


(75a) 

(75b) 

(76) 


with  Q (2yll)  laY  /dl  - 0.  Solution  of  this 
equation  is  sought  of  the  form 


x <*  .’J/'  £ Am,? 


•xp  - Z l£  iq 


m 0 


with  A 0 for  m odd  lor  the  is>  under  exan 
ination.m  Here,  for  simplicity,  only  the  m 0 
term  is  pursued  (with  A^'£),  hem  eforth,  lenot, 
by  A (£  ) ),  such  (hat 

X l(  ,t|'l  A (£  I exp  - / (£  iq  1 ■ . 

This  form  satisfies  (73b).  For  this  form  to 
satisfy  (73c),  it  is  determined  that 


On  the  left-hand  side  of  (73a),  the  approxi- 
mations tanh  q *>  q and  sech  q 1 are  employed: 
while,  on  the  right-hand  si^le,  the  approximation 
# I F*(q)|  Y0<£,q)  t»  aY0ooq  is  employed.  A con- 
ventional procedure  is  to  assign  the  constant  a by 
requiring  that  the  solution  to  the  substitute  equa- 
tion satisfy  the  original  equation  in  an  integral 
sense  over  the  range  of  independent  variables  of 
interest,  i.  e. , 1 s £ < 00  . 0 s q < °»  . This  pro- 

cedure is  unwieldy  in  the  present  case.  Here,  for 
explicitness,  it  is  taken  that  a may  be  assigned 
by  the  equation  of  aY^^q  to  £ | F I Yq  at  the 
value  of  q , for  which  the  latter  expression  is  a 
maximum.  (1)  If  it  is  arbitrarily  taken  that 
Yq~  Yo,.  then>  ilf'l  =<p(l  -«32)  . _with 
C - tanh  q , is  a maximum  at  q = tanh  * {1/^3  } , 
and  a # 0.888.  However,  such  a value  for  the 
factor  a is  too  large,  because  Yq  has  been  set 
to  its  maximum  value,  a value  attained  only  at 
large  £ and/or  q,  where  the  combustion  term 
is  sm?U.  (2)  if  it  is  taken  that  Yqw  Y 
Y 'l-(I-<0  )°  V (cf.  (51b)),  then,  a becomes 
a function  of  a , i.e.,  a = a(<7).  In  fa^t,  the 
equation  of  aq^  to  <p(  1 -<fl^)  j 1 - ( 1 -<fi  ) } at  the 
value  of  q for  which  this  latter  expression  is  a 
maximum  yields  a(£)  = 0.0928,  a(l>T  0.  387, 

...  . Such  values  for  a are  too  small,  because 
throughout  much,  though  not  all,  of  that  portion  of 
the  £,q-plane,  where  the  combustion  term  is 
appreciable,  Y_  > Yq  . However,  such  esti- 
mates suggest  trie  range  of  interest  for  values  of 
a for  physically  pertinent  values  of  the  other 
parameters.  Actually,  numerical  experimen- 
tation by  the  authors  indicates  that  a(&)  ± 0.28, 
a(l)  f 0,48,  ...  succeed  fairly  well.  However, 
the  adequacy  of  these  values  should  be  examined 
in  each  instance. 


A (0)  A = Y 

i h ) 


\ " 


l 


I 2 


* * 

Z (0)  = Z.  - -a  - j T i 
1 < V 2 ' (2(7 )*  H * 


( 78a  ' 


(78bl 


Sipce  H = H(<7),  Z.  =Z.(<7),  with  Z (J  ) 1 /ir  , 

Zi(l)=ir/8,  Z*(3/,2)  ^ 4/3rr,  Z { ( 2 ) l9*/64 

A (£  ) and  Z (£  )•  for  £ > 0 , are  assigned  to 

satisfy  the  approximate  partial  differential  equa- 
tion, (76). 

Substitution  of  (77bl  into  (76)  gives 

*'  * * 

[A  + Z A ] 

* *2  r *'  * 

-A  q [Z  (1  - 2Z  l-?Q  0,  (79) 

where  primes  denote  differentiation  with  respect  to 
{ . Equation  of  the  coefficients  of  powers  of  q to 
zero  yields 


* ' * # . :t 

Z - Z (1  - 2Z  ) = SQ  ; 

$ ' * * 

A + Z A = 0 . 


(80a  I 
(80b) 


The  solution  for  Z (£  ) is  determined  to  be 


* 

Z = 


(S*  4 1)  - (S  ‘ - 1)  F exp  ( - S £ J 


4 1 1 + E exp  { - S £ ] 


* * i 

S = (1  MQ  )z i 1 ; 

* * 

* (S  +U-4Z. 

=— 7 7 

(S  - 1)  +4Z. 


(81a) 


Since  (S  t 1)  ^ 2 and  since  4 /it  5 Z(  4 9*/  16  < 2 
for  i 4 <J  4 2 , it  follows  that  E > 0 (for  this 
range  of  rr ) . In  turn,  it  is  determined  that  the 
solution  for  A (£  ) is 


r4 


* * r r s * 

L = A.  exp  |-  \ Z • d£  J ' 


\ 


= A 


•>  i 

i [— . ~ ! "xpH,s 


I >E"  exp  f-S  £ i 
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Based  upon  these  solutions  for  Z and  A , it 
is  found  that  the  approximate  solution  for  X(4,  V ) 
takes  the  form 


XQ ' > 0 
XQ*  = 0 


V 


V 


> 1 
= I 


As  4 - 1 , 

X-  Y exp  |-$<T(4Z*)  T?2}  , (83a) 

with  Z.  given  in  (78b);  while,  as  £-  ® , 

, * 

X - Yf.  4*t(S  + " exp  j-J)(7  (S  " + l)r)2j  . (83b) 

From  inspection  of  this  approximate,  Gaussian- 
like  form  for  X,  it  is  noted  that  the  width  of  the 
profile  decreases  with  increasing  4>  for  fix^d  V. 
(This  result  holds  for  (S  + 1)  > 4Z* , i.  e. , E > 0.  ) 

* * 

For  Q =0  and/or  S =1,  i.  e. , for  no  chemi- 
cal reaction,  the  approximate  solution  of  (82) 
becomes 


For  this  case:  as  4 1 . 


X-  YFj  exp{-<j(2Z.  )rj2}  , (85a) 

just  as  for  the  general  case;  while,  as  £ -*®, 

. 1 7 

X-YFj4  exp{-<yr)  }.  (85b) 

The  self-similar  exact  solution  for  X for  the  case 
of  no  chemical  reaction,  as  calculated  from  (45) 
with  B = 0 , is 

X ««  Y 4 * (1  - tanh2r))<T  . (86) 


For  4-  ® , it  may  be  seen  that  the  approximate 
solution  decays  appropriately  in  f)  as  T)-*0,  but  too 
rapidly  as  TJ  -•  «,  Still,  the  general  behavior  is 
recovered. 

+ Division  of  (83b)  by  (85b)  indicates  that,  for 
Q > 0 and/or  S > 1 , i.e.,  for  finite  chemical 
reaction,  as  4 "•  * > 


«=  4 4 <S  " 1 ' exp  - ■ <T  (S  - 1 1 1)  ' . 1 87 1 

Thus,  for  finite  chemical  reaction,  the  fuel  pro- 
file tends  to  be  of  smaller  width  and  magnitude 
(than  for  no  chemical  reaction!. 

V 11.  Numerical  Solutions 


A conventional  C rank-Nicolson  implicit  finite- 
differencing  scheme  (cf.  Bush,  Feldman  and 
Fendell^)  is  employed  to  obtain  numerical  solutions 
to  the  boundary-value  problem  of  (451.  In  fact,  the 
program  has  been  written  to  solve  the  differential 
equation 


X,  + 2a<pX,  - 40  (1  -o"')4X,  , 

W V 4 


= (8cryB)  ((p(l  -o'"))  . 


Yq*"  c.4-r(.-o2»a-x  , 


r 


( 4 5a  l 


Clearly,  for  m = n = 1 , (45a) ' reduces  to  (45), 
i.e.,  this  model  with  the  more  generally  ordered 
reaction  expression  becomes  the  model  treated 
heretofore.  In  the  scheme  used  here,  three-point 
differences  approximate  the  derivatives  in  r; 

(for  0 * T)  s 9)  , with  the  step  size  prespecified 
and  constant  (i.e..  At)  = 0.  1);  while,  two-point 
differences  approximate  the  derivatives  in  4 
(cf.  Schlichting  7).  Extrapolation  of  rates  of 
change  in  4 of  previous  steps  establishes  the 
local  step  size  in  4 (except,  of  course,  for  the 
initial  step  in  4'*  The  even-symmetry  gradient- 
type  condition  at  the  axis  is  treated  (in  standard 
fashion)  by  the  simplest  central  difference  ex- 
pression. The  value  of  the  dependent  variable  at 
the  one  point  (external  to  the  field)  introduced  by 
this  procedure  is  eliminated  by  having  the  dif- 
ferential equation  satisfied  at  the  boundary  point 
(Fox^).  For  m =/  1 , Taylor  series  expansions 
are  used  to  linearize  the  nonlinear  (chemical  con- 
sumption) term  about  the  nominal  local  solution 
for  the  fuel  mass  fraction;  thus,  only  a simple 
inversion  of  a tridiagonal  matrix  is  required  for 
the  linear  algebraic  equations  arising  from  the 
finite-differencing.  Convergence  in  T)  dependence 
is  taken  to  be  effective  invariance  between  suc- 
cessive approximations. 


Through  numerical  integration,  using  (8a), 
(23a),  and  (50),  it  is  readily  possible  to  transform 
back  to  x,  y-coordinates  from  4.  ^-coordinates. 


It  is  convenient  to  characterize  results  in 
terms  of  the  integral  of  the  fuel  mass  fraction. 
Ip.  * Ij-(  4)  . and  the  normalized  integral  of  the 
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fuel  mass  fraction,  I ^ = I ^ (O  - )/l  j,(l ) = 

I_ / 1 , where 

F F i 

00  30 

I = [ Xdt);  T . = [ X.  dr).  (88) 

F F.  JQ  i 

Numerically,  1^.  is  obtained  by  application  of 
Simpson's  rule. 

VI  11.  Results 

A review  of  the  experimental  data  on  gaseous 
turbulent  fuel  jets  exhausting  into  air  is  useful, 
before  a presentation  of  the  predictions  of  the 
present  model  is  given.  In  this  way,  the  reader 
may  be  made  aware  of  what  hot-flow  measure- 
ments are  actually  available  to  validate  models. 

The  only  data  on  planar  turbulent  free-jet  dif- 
fusion flames  (known  to  the  authors)  are  the  very 
limited  measurements  due  to  Kremer,  who 
examined  mainly  axial  behavior  for  city  gas 
iprinc  ipally  hydrogen  and  methane)  issuing  from 
a 150  mm  slot.  In  100  slot  widths  downstream, 
the  axial  fuel  mole  fraction  dropped  by  approxi- 
mately one-half;  while,  the  axial  temperature 
increased  from  100°C  to  800°C.  A Reynolds 
number  of  3320,  based  on  slot  width,  produced 
a fully  turbulent  flame.  Also  qualitatively,  the 
spreading  of  the  reacting  flow  with  downstream 
distance  was  less  than  that  of  an  isothermal  jet. 
Kremer  recommends  use  of  a turbulent  Schmidt 
number  or  ==  0.75.  These  few  statements  are 
about  all  the  data  available  on  planar  jet  diffusion 
flames. 

Thus,  it  is  necessary  to  consider  experimental 
results  from  fuel  jets  issuing  from  round  orifices, 
even  though  one  can  only  qualitatively  compare 
axisymmetric  data  with  planar  theory.  It  may  be 
remarked  that  Bangert  and  Roach^  recently  com- 
pared predictions  (obtained  by  numerical  inte- 
gration), based  on  a parabolic  reactant-consump- 
tion model,  very  similar  to  the  one  proposed  here, 
to  the  data  (for  the  radial  profiles  of  temperature, 
fuel  mole  fraction,  oxidant  mole  fraction,  and 
product  mole  fraction)  furnished,  at  several  down- 
stream positions  for  an  axisymmetric  hydrogen 
fuel  jet  exhausting  into  coflowing  air,  by  Kent  and 
Bilger*^  (see,  also,  Bilger  and  Beck*®).  Bangert 
and  Roach  found  that,  for  m = n = 1 , B = 12  gave 
uniformly  excellent  agreement  between  prediction 
and  experiment  for  the  jet  speed/coflowing  ambient 
stream  speed  ratio  ranging  from  10  to  2;  further- 
more, there  was  rather  weak  sensitivity  to  mode r - 


Kremer's  conclusion  is  based  partly  on  non- 
reacting flow  data;  however,  for  axisymmetric 
turbulent  jet  diffusion  flames,  Elckhoff^'’  states: 
'Compared  to  the  data  evaluated  for  non-reacting 
jets  ...»  there  seems  to  be  no  significant  in- 
fluence of  combustion  on  the  ratio  of  turbulent 
momentum  to  mass  transport.  " 


ate  changes  in  B (say,  an  increase  nr  decrease 
in  magnitude  of  30%)  (Bangert,  private  communi- 
cation!. In  pa  rticula r , the  overlapping,  coexisting 
mean  profiles  for  the  reactant  mass  fractions,  well 
known  since  the  results  of  Hawthorne,  Weddell  and 
Hottel,  *^  are  recovered  by  the  model.  Also  re- 
covered by  the  model  are  (I)  the  ini  rease  and,  then, 
decrease  of  the  axial  temperature,  (2i  the  mono- 
tonic decrease  of  the  axial  fuel  mole  fraction,  and 
(3)  the  monotonic  increase  of  the  axial  oxidant 
mole  fraction,  all  with  increasing  downstream 
distance,  as  reported  by  Kent  and  Bilger.  Here, 
the  aim  is  to  point  out  other  experimental  evidence 
that  validates  the  model.  Furthermore,  the 
present  model  is  appreciably  simpler  than  that  of 
Bangert  and  Roach  in  that  ( 1 ) it  uses  a coordinate- 
transformation  to  correlate  the  compressible  and 
incompressible  equations:  (2)  it  uses  an  explicit 
eddy  viscosity  in  preference  to  a two-equation 
field-type  closure  to  specify  the  turbulent  diffusion; 
and  (3)  it  uses  a facile  self-similar  flow  field  (for 
a stagnant  ambient)  in  preference  to  calculating 
the  actual  nonself-similar  flow  field.  These 
appreciable  simplifications,  previously  shown  to 
permit  approximate  closed-form  analytic  charac- 
terization of  the  model,  are  revealed  below  to 
sacrifice  little  in  the  way  of  predictive  capacity. 

For  axisymmetric  turbulent  fuel  jets  exhausting 
into  stagnant  air  (for  cases  of  high  jet  speed  and 
small  diameter,  such  that:  the  Reynolds  number  is 
large,  so  that  molecular  processes  are  negligible 
relative  to  inertial  effects;  and  the  Froude  number 
is  large,  so  that  buoyant  accelerations  are  negli- 
gible relative  to  inertial  accelerations),  the  fol- 
lowing experimental  results  are  reported. 

1.  Reynolds  number.  Not  only  is  the  length  of  the 
turbulent  flame  known  to  be  but  weakly  a function 
of  the  Reynolds  number  (Hawthorne,  Weddell  and 
Hottel*"’),  but  also  the  axial  temperature  distri- 
bution is  but  weakly  dependent  on  the  Reynolds 
number.  Kremer*^1  states  that  the  number  of 
nozzle  diameters  downstream  at  which  the  maxi- 
mum temperature  occurs  is  increased  by  25%,  and 
that  the  magnitude  of  the  maximum  axial  tempera- 
ture is  increased  by  15%,  when  the  Reynolds 
number  is  increased  by  a factor  of  4. 

2.  Temperature  distribution.  At  fixed  small 
distances  downstream,  the  maximum  temperature 
occurs  several  nozzle  diameters  off  the  axis 
(Chigier  and  Strokin^).  The  radial  position  of  the 
maximum,  and  the  magnitude  of  the  maximum, 
first,  increase,  then,  decrease  with  increasing 
axial  distance  downstream.  The  magnitude  of  this 
maximum,  with  respect  to  downstream  position, 
rises  but  slowly  after  about  10  nozzle  diameters 
downstream;  peaks  at  a value  below  the  adiabatic 
flame  temperature  at  about  that  downstream  dis- 
tance at  which  the  maximum  goes  to  the  axis 
(roughly  100  diameters);  and,  then,  decreases 
appreciably  with  farther  downstream  distance 
(Kremer*^;  Lockwood  and  Odidi5®).  I he  tempera- 
ture on  the  axis  rises  monotonically  until  about 
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100  diameters  downstream;  then,  decreases 
IKremer^;  Guenther  and  Simon*’).  These  re- 
marks concerning  the  temperature  field  hold  for 

a coflowing  oxidant  stream  as  well,  whether  that 

40 

stream  be  heated  llaleno  and  Kotani  ) or  not 
(Kent  and  Bilger*7;  Bilger  and  Beck'®). 

3.  Flame  contour  [locus  of  equal  (stoichiometri- 
cally  adjusted)  mass  fractions  for  fuel  and  oxidant, 
where,  experimentally,  an  effective  composite 
fuel  profile  is  defined].  As  interpreted  by  Zimont 
and  Meshcheryakov,  4 Kremer**'  finds  that  the 
flame  contour:  first,  increases  with  downstream 
distance;  then,  decreases  to  the  axis.  However, 
the  flame  contour  envelopes  the  maximum  tem- 
perature contour  (i.e.,  the  flame  contour  lies  at 
a larger  radial  position  at  all  downstream  stations, 
and  goes  to  the  axis  farther  downstream  than  does 
the  maximum  temperature  locus).  Measurements 
of  the  distance  downstream  at  which  the  flame 
contour  goes  to  the  axis,  i.e.,  the  flame  length, 
are  given  by  Kremer,"’who  found  that  for  other 
parameters  held  nearly  constant,  a city-gas  (H^, 
CH|,  CO,  CnHm;  N'2.  CO2)  turbulent  jet  dif- 
fusion flame  is  one-and-one-half  times  as  long  as 
a carbon  monoxide  flame,  and  a hydrogen  flame  is 
two-and-one-quarter  times  as  long.  Most  other 
purported  measurements  are  based  on  visual 
observations  (cf.  , e.  g.  , Hawthorne,  Weddell  and 
Hottel**)  and  are  not  reported  here.  Incidentally, 
Brzustowski4^  finds  that  published  analytic  ex- 
pressions, asserted  by  their  proponents  to  be 
adequate  lor  predicting  flame  length,  are  not 
adequately  verified.  This  is  not  unexpected,  in 
that  such  formulas  include  unattained  quantities, 
such  as  the  adiabatic  flame  temperature  (e.g., 
Hawthorne,  Weddell  and  Hottel**),  or  ill-defined 
quantities,  such  the  mean  flame  density  (e.  g.  , 
Beer  and  Chigier  ). 


lft 


4.  Species  profiles.  On  the  axis,  the  fuel  mass 
fraction  monotonically  decreases,  while  the  oxi- 

dant mass  fraction  monotonically  increases,  with 
increasing  distance  downstream.  The  product 
mass  fractions  increase  until  the  fuel  is  virtually 
exhausted,  and  the  maximum  axial  temperature  is 
achieved;  then,  the  product  mass  fractions  de- 
crease with  further  downstream  distance  (Kremer 
Guenther  and  Simon*^).  At  a fixed  downstream 
station,  the  fuel  mass  fraction  decreases  mono- 
tonically and  the  oxidant  mass  fraction  increases 
monotonically,  with  increasing  radial  distance 
from  the  axis.  At  a shorter  downstream  distance, 
the  product  mass  fractions  increase  to  a maximum 
a few  nozzle  diameters  off  the  axis;  then,  decrease. 
At  about  the  same  downstream  distance  at  which 
the  fuel  species  are  effectively  totally  depleted  and 
the  maximum  temperature  is  first  achieved  on  the 
axis,  the  maximum  of  the  product  mass  fractions 
is  also  achieved  on  the  axis  (Hawthorne,  Weddell 
and  Hottel*®:  Kremer  Lavoie  and  Schlader4*). 
These  same  characterizations  hold  for  a coflowing 
ambient  oxidant  stream  as  well  (Kent  and  Bilger*  , 
Bilger  and  Beck*®).  Chigier  and  Strokin'*7  show 
that  a passive  scalar,  such  as  the  mass  fraction 


for  an  inert  species,  in  an  axisymmetric  turbulent 
jet  diffusion  flame  employing  methane  in  stagnant 
air  at  an  exit  Reynolds  number  of  hnOO,  satisfies 
similarity  concepts,  in  that  the  passive-scalar 
field  ( I I decays  (mono ton ica  1 1 y I axia  11  y . w it h t hi  s 
decay  inversely  proportional  to  the  first  power  of 
the  axial  distance;  (21  decays  I monotmi ica  1 1 y ) 
radially,  with  this  decay  a (Gaussian1  function  of 
the  square  of  the  radial  distance/axial  distance 
ratio.  The  rate  of  axial  decay  is  about  0.4  1 times 
as  rapid  as  that  for  the  nonreacting  case. 


5.  Rate  of  spread.  Turbulent  axisymmetric  iet 
diffusion  flames  have  long  been  known  to  have  a 
significantly  lower  rate  of  spread  than  a constant- 
density  nonburning  jet  (Beer  and  C hig le r “ * ; 

Chigier  and  Strokin5'!.  Whereas  a nonburning  iet 
may  spread  at  a half-angle  of  5-7.  5°,  a burning 
jet  spreads  typically  at  2.  5-3°.  As  noted  by 
Zimont  and  Meshcheryakov,4*  confirmatory  evi- 
dence is  presented  by  Kremer.  * In  fact,  in  a 
related  statement,  Zimont  and  Meshcheryakov 
note  that,  in  a nonreacting  jet,  the  analogous  con- 
tour for  the  analogous  passive  scalar  (linear 
coupling,  or  Shvab-Zeldovich,  functionl  is  uni- 
formly enveloped  by  the  flame  contour,  and  the 
axial  intercept  in  the  nonreacting  case  is  very 
nearly  half  the  distance  downstream  that  holds  for 
the  reacting  case.  Effectively,  then,  > is  reduced 
for  the  reacting  case.  Equivalently,  x^  is  re- 
duced, since  x;  and  Z are  linearly  proportional 
through  z^,  i.e 
modifying  Xj  , 
empirical  adjustment,  as  opposed  to  modifying  y , 
in  a rate-of-spread- t y pe  empirical  adjustment, 
is  set  forth,  for  example,  by  Vulis  and  Yarin.'*'*) 


. , Xj  = z.y  . (The  concept  of 
in  a contractual-mapping-type 


6.  Velocity  Field.  Chjgier  andStrokin  'and  Zimont 
and  Meshcheryakov  furnish  particularly  clear 
comparisons  between  nonreacting  and  reacting 
jets,  for  axial  decay  of  the  velocity  and  dynamic 
pressure.  The  decay  from  initial  value  is  re- 
duced in  the  reacting  case. 


7.  Turbulent  diffusion  coefficient.  Whereas  the 
previously  discussed  quantities  may  be  directly 
measured,  ||^is  quantity  must  be  deduced.  C Inkier 
and  Strokin  and  Zimont  and  Meshcheryakov 
concur  that,  axially,  the  turbulent  diffusion  coef- 
ficient for  the  reacting  case  is  smaller  than  that 
for  the  nonreacting  case  for  about  10  nozzle 
diameters  downstream;  however,  by  50  nozzle 
diameters  downstream,  the  reacting-jet  coefficient 
is  twice  as  large,  and,  by  100  nozzle  diameters, 
three  times  as  large.  Chigier  and  Strokin*'  also 
suggest  that  the  turbulent  diffusion  coefficient 
decays  to  one-third  its  axial  value  at  large  radial 
distances,  at  distances  of  20-30  nozzle  diameters 
downstream. 


Attention  is  now  turned  from  exposition  of 
experimental  data  to  presentation  of  predictions  of 
the  model.  First,  however,  it  is  noted  that,  em- 
pirically, y 7.67  for  the  planar  nonreacting  case 
(Schlichting  ’);  it  is  known  that  y is  reduced  in 
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value  for  the  chemically  reacting  case  (Kremer  ). 
Here,  based  upon  axisy mmet r ic  data,  since 
quantitative  plana  r data  are  not  available,  7 is 
set  equal  to  one-half  of  its  inert  value.  Second, 
the  specific  case  chosen  for  study  is  a planar 
turbulent  jet  of  pure  carbon  monoxide  gas  at  room 
temperature  exhausting  into  stagnant  air,  with 
the  reaction  taken  to  proceed  by  a direct  one-step 
irreversible  mechanism  to  form  the  product  gas 
ca  rbon  dioxide.  Accordingly,  parametric  values 
are  assigned  as  follows  : Y_  = 0.  638,  Y - 
1.56,  YQ,  0,  T.  = 1,  Q = and*--  3^35. 

The  parameter  cr1  is  to  be  set  equal  to  either  £ 
or  1;  while,  yB  is  to  be  set  equal  to  0,  0.  1,  or 
I,  for  purposes  of  comparison.  Finally,  in  the 
expression  for  mean  reactant  consumption,  the 
exponent  on  the  fuel  mass  fraction  m is  to  be 
set  equal  to  either  1 or  2,  for  parametric  study; 
while,  the  exponent  on  the  oxidant  mass  fraction 
n is  to  be  set  equal  to  1 . Clearly,  the  vari- 
ability in  assignment  of  values  for  <7  , yB,  and 
m reflects  uncertainty  concerning  appropriate 
values  of  these  parameters. 

Some  results  may  be  noted  prior  to  the  compu- 
tations. First,  for  a pure  fuel  jet  exhausting  into 
invitiated  air,  the  flame  length  4f  jl  + (128l£j)/ 
(m^  )}  . For  hydrogen,  Uq=  1 , = 2 , 

m^,  = 2;  while,  for  carbon  monoxide,  Pq  = 1, 

P^.  = 2,  m^.  = 18.  Thus,  the  flame  length  for 
hydrogen  is  appreciably  longer  than  that  for  carbon 
monoxide,  although  quantitative  comparison  is 
limited  by  the  fact  that  the  approximation  of  com- 
parable molecular  weights  for  all  species  has  been 
adopted.  In  any  case,  it  may  be  premature  to 
ascribe  difficulties  in  predicting  turbulent  jet 
diffusion  flame  lengths  to  the  neglected  role  of 
molecular  diffusion  (Spalding13),  as  opposed  to 
inadequate  modeling  of  large-scale  mixing.  Also, 
since  4 = x/xj  = (4Tj/3)  (x/y)  , it  becomes  evident 
that  the  empirical  decrease  in  y implies,  within 
the  context  of  the  present  model,  a decrease  in 
the  axial  decay  of  the  velocity  and  of  a passive 
scalar  (cf.  (32b)  and  (41),  respectively).  Further, 
since  < = DT2  = (3/l6)(4T:/3)*(x/y3)*T2  , the 
empirical  decrease  in  y,  but  eventual  increase 
in  T owing  to  chemical  exothermic  ity,  implies, 
within  the  context  of  the  present  model,  that  the 
turbulent  diffusion  coefficient  for  reacting  flows  is 
modestly  smaller  than  for  frozen  flow  at  small 
downstream  distances,  but  appreciably  larger  for 
large  downstream  distances.  Thus,  the  present 
model,  for  empirically  assigned  y,  can  recover 
several  known  trends  in  the  data. 

Consideration  of  the  frozen  chemistry  case, 
with  B - 0,  permits  comparison  of  the  numerical. 


It  is  recalled  that  Y = (m  Y *)/(m  V ) , 

U U U O 

Y - (mY  )/(m  V ),  where  m = m P + 
F f + F F O O 

m^Pj,  , and  Y^  is  the  mass  fraction  of 

species  1. 


iterative,  and  Oseen-linearized  solutions  with  the 
known  exact  solutions,  for  error  analysis.  For 
B O,  the  mass  fractions  for  fuel  and  oxidant, 
individually,  are  passive  scalars,  with  constant 
flux  integrals  (cf.  (36)  ).  1-  or  B 0 , it  is  ex- 
pected that  Njl?!  N,(4>.  for  4 1.  where 

CD 

N 1 = \ (1  - tanh^tj)  • X(4,Tj)dtJ: 

' 0 

3D 

- 5 ? 

N (41  - 4 (l-tanh  'll1  Xll.ni  dO  . 189! 

2 0 

The  accuracy  of  the  numerical  integration  is 
characterized  by  the  results  that,  for  a 1 , at 
4 = 15,  N = 0.2696,  N 0.  2b78;  while,  at 
4 = 140,  N 0.08886,  N 0.08764.  The  axial 
values  for  the  fuel  mass  fraction,  X(4,0l=  X°  4), 
found  by  the  iteration  scheme  are  uniformly  less 
than  1%  different  from  (smaller  than)  those  found 
by  the  numerical  integration,  for  1 * 4 4 190. 

The  values  for  X°(4*  found  by  the  Oseen-linear- 
ization  method  range  between  4 ' and  1(6  below 
those  found  by  the  numerical  integration,  again, 
for  1 s 4 s 140.  For  <7  ?,  the  flux  integrals 

incur  about  the  same  percentage  discrepancy  as 
for  <7  = 1;  however,  the  axial  values  found  by 
iteration  are  uniformly  less  than  117  too  low,  while 
the  Oseen-linea  rization  values  are  up  to  2017  too 
high,  relative  to  X°(4)  obtained  numerically,  for 
I s 4 s 140,  (Clearly,  adjusting  the  Oseen  constant 
a does  not  avail  for  B = 0,  and,  in  fact,  intro- 
duction of  a second  Oseen  constant  with  respect  to 
analytically  motivated  modifications  of  the  flow 
field  might  have  been  undertaken  for  purposes  of 
better  recovery  of  the  fuel  mass  fraction  field  in 
the  frozen  case;  however,  this  procedure  was  not 
adopted.  ) 

Results  for  the  nominal  case,  i.e,,  yB  1, 
m = 1,  n = I,  <7  = $,  are  plotted  in  Figs.  2-11. 

For  rj  = 0 , the  monotonic  decrease  of  Y (and 
PU2)  , and  the  monotonic  increase  of  Y^,  with 
increasing  4.  are  depicted  in  Fig.  2.  The  mono- 
tonic decrease  of  Y , and  the  monotonic  increase 
of  Yq,  with  increasing  t? , at  any  fixed  4 * 1 . 
are  depicted  in  Fig.  3.  The  maximum  temperature 
occurs  off  the  axis  for  4<4f  11.82  (the  flame 

tip  distance),  but,  then,  goes  to  the  axis  for  larger 
4 . As  seen  in  Fig.  4,  the  magnitude  of  the  maxi- 
mum temperature  increases  with  increasing  4 for 
4 < 4°  , though  most  of  the  increase  occurs  at  rela- 
tively small  4 ! then,  this  magnitude  decreases  for 
4 > Cj  • These  results  are  amplified  in  Tig.  5, 
where  it  is  shown  that,  at  given  4,  the  y-position 
of  the  locus  Yp  = Yj.  = increases,  then,  de- 
creases, with  increasing  4<  hut,  nevertheless, 
always  exceeds  the  y-position  of  the  maximum 
temperature.  All  these  results  are  compatible  in 
trend  with  experimental  data,  discussed  earlier, 
and  are  of  plausible  magnitude,  when  allowance  is 
made  for  the  planar  geometry.  Whereas  Figs.  2-6 
involve  the  numerical  solutions  only.  Figs.  6-10 
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indicate  the  adequacy  of  the  approximate  methods 
of  solution.  (Mere,  for  O f,  the  Oseen-l  inea  r- 
ization  constant  a is  set  equal  to  0.  38.  as  opposed 
to  the  0.28  value  suggested  earlier  after  numerical 
experimentation.)  Fig.  11  reveals  that,  in  the 
conservation  equation  for  fuel,  turbulent  diffusion 
and  tangential  (streamwisel  convection  balance 
near  the  axis,  but  normal  convection  and  chemical 
reaction  (consumption)  play  a significant  role  off 
the  axis;  at  large  distances  from  the  axis,  diffusion 
and  normal  convection  dominate.  It  seems  worth 
noting  that  plausible  profiles  are  generated  by  a 
model  that  uniformly  permits  no  consumption  on 
the  axis. 

In  Table  1,  numerical  results  for  the  nominal 
case  are  compared  with  numerical  results  obtained 
from  the  model  by  sequentially  altering  one  param- 
eter at  a time.  Increasing  the  turbulent  Schmidt 
number  increases  the  axial  values  of  the  fuel  mass 
fraction  and  the  distance  off  axis  of  the  flame. 
Increasing  the  fuel  or  oxidant  exponent  in  the  con- 
sumption term  decreases  the  amount  of  fuel  con- 
sumption: while,  increasing  the  effective  turbulent 
Damkohler  number  yB  increases  the  amount  of 
fuel  consumption.  (By  numerical  experimentation 
for  the  case  of  yB  = (7=m  = n=l,  in  the  Oseen- 
linearization  model,  a = 0.62,  as  opposed  to 
a = 0.48,  mentioned  above,  gives  closest  agree- 
ment with  the  results  from  the  finite-difference 
integration;  in  fact,  the  axial  fuel  mass  fraction  is 
missed,  i.  e.  , over  estimated,  by  less  than  1%, 
for  1 * i ‘ 30.  ) 

In  conclusion,  the  authors  feel  that,  in  the  pres- 
ent paper,  the  following  have  been  demonstrated: 


(1)  the  postulated  model  captures  the  essential 
physics  (i.e.,  the  experimental  results!  for  the 
turbulent  portion  of  the  planar  fuel  let;  and  (2!  the 
given  approximate  solutions  capture  (lie  esscnt  e of 
the  model.  It  is  suggested  that  studies  of  other 
parametric  cases  may  be  made  using  the  approxi- 
mate solutions,  rather  than  the  numerical  solutions. 
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Fig.  1.  Schematic  diagram,  not  to 
scale,  of  a lifted  subsonic 
planar  turbulent  fuel  jet 
exhausting  into  a stagnant 
oxidant-containing  ambient, 
with  the  flame  locus  defined 
by  the  condition  of  equal 
values  for  the  (stoichiometri- 
cally  adjusted)  mass  fraction 
for  fuel,  Yr,  and  for  oxidant, 
Yq.  At  the  axis  of  symmetry 

y = n = 0;  x,  £ are  streamwise 
coordinates. 


Fig.  2.  For  the  nominal  case  (Yq  = 

oo 

0.638,  Yc  = 1.56,  Y*  = 0, 
hj  Ui 

T,  - 1,  Q = 64.4,  y = 3.835, 

By  = 1 , o = 0.5,  m = n = 1 ), 
the  following  quantities  are 
plotted  against  the  stream- 
wise  similarity  coordinate 
t,  from  the  break  point 
£ * 1:  the  axial  fuel  mass 

fraction,  YpU.O);  the  axial 

oxidant  mass  fraction, YQU, 0) ; 

the  Integral  of  Yp(s,n)  over 

the  transverse  similarity 
coordinate  n,  IpU);  and  the 
axial  dynamic  pressure,  pU,0)* 
U2U.0). 


Fig.  3.  For  the  nominal  case,  the 

static  temperature  T ( • ) is 
qiven  as  a function  of  the 
(density-scaled)  transverse 
similarity  coordinate  n,  at 
several  values  of  the  stream- 
wise  similarity  coordinate  , 
where  £ = = 11.82  is  the 

flame  tip  [Yp(cJ  ,0)=YQ(r°  ,0)]. 


Fig.  4.  For  the  nominal  case,  the  fuel 
mass  fraction  Yp  (solid  curves) 

and  the  oxidant  mass  fraction 
Yq  (dashed  curves)  are  given  as 

a function  of  the  transverse 
similarity  coordinate  n»  at 
several  values  of  the  stream- 
wise  similarity  coordinate 
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Fig.  5.  For  the  nominal  case,  the 

following  quantities  are  given 
as  a function  of  the  streamwise 
similarity  coordinate  5,  for 
i 5 1 (where  £ = 1 is  the  break 
point):  the  value  of  Yq(=  Yp) 

at  the  flame  locus  for  £ < 11.82 
(where  5 = 11.82  is  the  flame 
tip  position),  this  curve  giving 
the  axial  value  of  the  oxidant 
mass  fraction  YqU.O)  for 

i > 11.82;  the  value  of  n at 
which  Yq  = Yp  for  given  5;  the 

equivalent  value  of  y at  which 
Yq  = Yp  for  given  £;  the  maxi- 
mum value  of  the  static  tem- 
perature T at  any  n for  given 
£ , Tj^;  the  value  of  n at 

which  occurs  for  given  £; 

and  the  equivalent  value  of  y 
at  which  occurs  for  given  £. 
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Fig.  6.  For  the  nominal  case,  a compari- 
son of  the  fuel  mass  fraction 
profiles  Yp  over  transverse  sim- 
ilarity coordinate  n,  for  £ = 1 
(the  break  point),  as  obtained 
by  numerical  integration,  by 
iteration,  and  by  Oseen  lin- 
earization (a  = 0.38). 


Fig.  7.  Analogous  plot  to  Fig.  6,  for 
C = 3.0. 
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Fig.  10. 


Analogous  plot  to  Fig.  6,  for 
£ = 15.0. 


Fig.  8.  Analogous  plot  to  Fig.  6,  for 
5 = 8.0. 


n 


Fig.  9.  Analogous  plot  to  Fig.  6,  for 
£ = 11.82  (the  flame  tip 
position). 


Fig.  11.  For  the  nominal  case,  at  the 
distance  downstream  of  the 
flame  tip  £;  = 11.82,  the  values 
of  the  terms  comprising  the 
model  mean  equation  of  conser- 
vation of  the  fuel  species 
[tangential  (streamwise)  con- 
vection, transverse  convection, 
turbulent  diffusion,  and  con- 
sumption (chemical  reaction)] 
are  given  as  a function  of  the 
transverse  similarity 
coordinate  n. 
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